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Abstract. The spreading of a circular liquid drop on a solid substrate can be 
described by the time evolution of its base radius R{t). In complete wetting 
the quasistationary regime (far away from initial and final transients) typically 
obeys the so-called Tanner law, with R ~ t"T , ax = 1/10. Late-time spreading 
may differ significantly from the Tanner law: in some cases the drop does not 
thin down to a molecular film and instead reaches an equilibrium pancake-like 
shape; in other situations, as revealed by recent experiments with spontaneously 
spreading nematic crystals, the growth of the base radius accelerates after the 
Tanner stage. Here we demonstrate that these two seemingly confiicting trends 
can be reconciled within a suitably revisited energy balance approach, by taking 
into account the line tension contribution to the driving force of spreading: a 
positive line tension is responsible for the formation of pancake-like structures, 
whereas a negative line tension tends to lengthen the contact line and induces 
an accelerated spreading (a transition to a faster power law for R{t) than in the 
Tanner stage). 
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1. Introduction 

The spreading of a liquid on a solid surface is a complicated process where many 
factors come into play, not necessarily known and not always controllable: the kinetic 
behavior may be strongly influenced or even dominated by the volatility and viscosity 
of the liquid (or by other rheological parameters if the liquid is non-Newtonian), 
by the presence of impurities in the bulk phases (chemical contaminants, surfactants, 
polymers, etc), by the roughness or texture of the surface, or by its crystalline structure 
and chemical composition. In consequence, the spreading of thin films is generally 
dependent upon the details of the structures and interactions in the co-existing phases. 
By contrast, for thicker films and drops one expects a less specific behavior, described 
by universal laws. However, although the prominent features of the spreading of 
macroscopic drops are relatively well understood [H HI [3] , a comprehensive theoretical 
framework in which all of the experimentally observed phenomena harmoniously find 
their place is still lacking at present. 

Here we are concerned with a standard textbook problem - the spontaneous 
spreading of a non-volatile drop on an ideal, fiat, clean, horizontal, homogeneous 
solid surface (see figure 1). The drop has a macroscopic size but is sufficiently 
small that we can completely discard the effects of gravity. The spreading parameter 
S — CTsG — 0-SL — cr is positivc, so that the drop tends to cover as much of the solid 
surface as possible to shield it against the gas phase; asG, ^sl and a are the interfacial 
tensions of the solid/gas, solid/liquid and liquid/gas interfaces, respectively. 

We focus on the following three known features of spontaneous spreading: 

• The base radius R{t) of a circular drop grows during spreading and, at 
intermediate times t (far from initial and final transients) , typically obeys a power 
law R{t) ~ t"^ with ax = 1/10 [4l[5], known as the Tanner law. 

• The final stage of spreading of non-volatile droplets is not always a molecular 
film. Sometimes a fiat, bounded structure is reached instead - a so-called pancake 
(see [T] and [eliaiH]). 

• An acceleration of the spreading process (an apparent transition from Tanner's 
power law to a faster one) has been observed for spontaneously spreading nematic 
liquid crystals [51 [lO]. Experiments revealed an algebraic growth R{t) ^ with 
a nearly twice as large as the exponent ax characterizing the Tanner law: a = 0.2 
[9] and a = 0.19 [10]. 

At present time, it is well understood why the radius R{t) of a spontaneously spreading 
circular drop grows in proportion to i^/^". This law has been derived analytically 
[H m [11] and confirmed experimentally on many accounts [U [HI HS] [14]. The 
fundamental argument is that the hydrodynamics in the bulk of a drop are driven 
by capillary forces alone, which directly yields R ~ for a self-similar bulk, in 

the lubrication approximation [H \5\. Alternatively, the trend can be regarded as 
a competition between the hydrodynamic dissipation (primarily in the contact line 
region of the drop) and an unbalanced capillary force [U [2 [3l [11]. Note that the 
law is rather universal - in the sense that the observed exponent 1/10 is most often 
independent of the precise nature of the spreading liquid - and has been observed 
not only for simple liquids, but for oils, polymeric liquids, liquid metals and nematic 
liquid crystals. For non-Newtonian liquids some deviations from the Tanner law can 
be observed, attributable to their specific rheological properties and hence, specific 
features of the hydrodynamic dissipation in the bulk. 
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Figure 1. a) A sketch of a macroscopic liquid droplet spreading on a solid 
substrate. R{t) and 9{t) denote the base radius and the contact angle of the 
macroscopic part of the droplet, respectively, as inferred from the inflection point 
at the apparent contact line, b) A sketch of a liquid droplet at a late spreading 
time, in the situation where the flnal stage of spreading is a mesoscopic pancake. 
An inflexion point may still exist, but the relevant R{t) now corresponds to the 
limit between the microscopic and mesoscopic regions; the corresponding 0{t) is 
undefined. 



The reason why in some cases a spontaneously spreading droplet attains an 
equilibrium pancake-like form is also clear. Such flat pancake-like structures are 
sometimes more favorable energetically than molecular films: this occurs when short- 
range interactions promote dewetting, even though the overall situation is that of 
complete wetting. Theoretically, such structures have been predicted and analyzed in 
p-J and [GlITj- They were also observed experimentally (see, e.g., [Hj). A key feature 
of the spreading process as it reaches a pancake shape - as a transient following the 
Tanner stage - is that the base radius of the drop tends to a stationary value (see figure 
lb). It must be noted, though, that the shape of such a drop differs considerably from 
that of a capillary cap; the geometrical meaning of R{t) is also quite different. 

As opposed to the Tanner law and the emergence of mesoscopic pancakes, the 
physical origin of the accelerated spreading observed in [9l [TOj has yet to be clarified. 
Also, the latter trend is apparently in conflict with the notion that a Tanner stage 
must be terminated by the onset of either a molecular fllm or a mesoscopic pancake. 

Our motivation in this paper is to identify the essential factors of the spreading 
kinetics that might have been disregarded so far, and thus to achieve a more 
complete qualitative understanding of the standard textbook problem. We focus more 
speciflcally on the spreading dynamics that may develop after the Tanner stage and 
try to account for the two apparently conflicting trends, possibly by introducing new 
notions or mechanisms. As a flrst step in this direction, we perform an analysis within 
the classical framework of an energy balance approach [1]. This approach is known 
to provide a qualitative derivation of the Tanner law, capturing only a few essential 
features of the phenomenon, in a physically transparent fashion. As noted in [T6|. 
the energy-based equations are functionally equivalent to the standard hydrodynamic 
approaches used in the literature but are lighter in terms of analytical calculations and 
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assumptions involved (as compared to, e.g., phenomenological boundary conditions, 
should the thickness profile be described by a differential equation). 

Within this approach, the Tanner law is obtained by balancing the rate of energy 
dissipation in the spreading macroscopic droplet and the driving force of spreading, 
which is taken equal to the unbalanced Young force. We point out that the analysis of 
de Gennes [T] disregards the Hne tension contribution to the driving force of spreading 
[l71[T8]. Typically, the line tension r is very small - only 10~^° to 10~^^ N (see, e.g., 
[191 l20l [21] ) - and it is legitimate to neglect it when deahng with large, essentially 
capillary droplets. Nonetheless we show that, when such a contribution is taken into 
account, a consistent, non-conflicting picture emerges, with the following trends for a 
drop (supposedly well-approximated by a spherical cap): 

• At sufficiently early spreading times the effect of r is negligible and the Tanner 
stage holds. 

• At long times and for negative values of r, the spreading process crosses over 
to a significantly faster power law than Tanner's R ^ t^/^° (as observed for the 
nematic droplets [9lfT0|). 

• At long times and for positive values of r, the growth of R{t) slows down and 
terminates at a finite value i?(oo). This latter trend is indicative of the emergence 
of pancakes. 

Therefore, the approach presented in this paper resolves a seemingly controversial 
behavior of spreading processes, provided that a properly defined line tension r is 
taken into account. We note, however, that this approach is justified only in the case 
of macroscopic drops, for which both the surface tension and the line tension are valid 
notions. 

The paper is outlined as follows. In section 2 we first present the derivation of 
the Tanner law in terms of the energy balance approach. Then, in section 3, we revisit 
the standard picture by analyzing different factors which may infiuence the spreading 
kinetics, especially the notion of line tension. Finally, in section 4 we discuss possible 
limitations of our approach. 

2. Energy balance approach 

To lay the basis of our analysis, we start with the derivation of the Tanner law within 
the framework of the energy balance approach, originally presented by de Gennes 
in his 1985 review paper [T]. Additional details, discussions and appHcations of this 
approach can be found in |14[ I16j. Further on, in section 3, we will discuss a couple 
of additional factors which are missing from the seminal approach, but may account 
for the abnormal spreading behavior observed after the Tanner stage. 

In figure la we sketched a typical configuration for a macroscopic Hquid droplet 
spreading spontaneously on a solid substrate. The drop can be "divided" into the 
following three regions: a "macroscopic" bulk, a "mesoscopic" film (a region that is 
within the range of surface forces), and a "microscopic" precursor, the thickness of 
which amounts to several molecular diameters. Note that figure 1 is schematic and 
the relative sizes of these three regions are not up to scale. For brevity, we will 
henceforth not use the adjectives "macroscopic", "mesoscopic" and "microscopic"; we 
will instead refer to the three regions as the bulk droplet, the film and the precursor. 

We note first that the edge of the precursor spreads well ahead of the film and 
the bulk droplet: the radius of the precursor grows as Vi [22 } [23 t [24 t [25 l [26 l [27 l [28] 
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and the film plays the role of a reservoir feeding the precursor; this picture is valid as 
long as the said reservoir is far from being exhausted. Thus the precursor is decoupled 
from the rest of the drop and its only role in the process can be seen as lubricating 
the substrate for spreading of the film and the bulk drop. 

A fundamental assumption of this approach is that the length scales of the bulk 
and the film are well separated, i.e., that the bulk of the drop is much wider and taller 
than the film. Thus the bulk can be adequately approximated by a spherical cap with 
the base radius R{t), contact angle 9{t), and nearly constant volume V {V = 
for sufficiently small 9), i.e., the bulk is in equilibrium at constant volume V and an 
instantaneous base radius R{t). 

Then, we can define an instantaneous free energy 2TrR{t)F{t) along with an 
instantaneous rate of energy dissipation 2TTR{t)W{t), and propose that the evolution 
of these two quantities obeys the standard relationship of the mechanics of dissipative 
systems: 

(W^macro + l^^mcso + l^micro) R = -U^^^^, (1) 

where U — dR/dt is the instantaneous velocity of the apparent contact line. The 
right-hand-side (r.h.s.) of HI) is the rate of change of the free energy of the system. It 
is equivalent to the power of the driving force dF/dR applied to the moving contact 
line, and is balanced by the total dissipation that occurs in the system, i.e., the left- 
hand-side (l.h.s.) of II]). The terms Wmacro, Wmcso and Wmicro are the dissipation 
rates in the bulk, film and precursor, respectively, divided by the length 27ri? of the 
apparent contact line. Keep in mind that F{t), W{t) and the driving force dF/dR are 
also, by definition, reduced by the length 2TrR{t). 
Next we specify the r.h.s. of ^ [l]: 

Equation ^ takes into account the surface energies of the three interfaces meeting at 
the macroscopic contact line and determines their variation with respect to R(t) at 
constant spherical cap volume. Note that the result is equivalent to a straightforward 
application of the Young law; in fact, the r.h.s. of ^ is typically referred to as an 
"unbalanced Young force". 

As for the l.h.s. of |(2l), we can formally decompose the dissipation according to 
the regions outHned in figure [IK, as follows: Wmacro corresponds to hydrodynamic 
dissipation in the bulk drop, where viscous fiows are driven by the capillary pressure; 
Wmeso corresponds to hydrodynamic dissipation in the film, where viscous fiows are 
driven by the disjoining pressure; Wmicro corresponds to friction at the microscopic 
scale, both at the edge of the film and in the molecular precursor. 

The dissipation in the bulk drop is well-approximated by that in a wedge, and is 
of the form 



Wmacro^ 3 ?7C/2^(0) In 



(3) 



where rj is the viscosity, Xmax and Xmin are effective cutoff lengths for the integration 
over the droplet height and g{9) is a known function of the instantaneous contact 
angle. A salient feature is that the leading asymptotic behavior of g{9) when 6* — > 
is g{9) ~ 1/9, which means that Wmacro exhibits an unbounded growth as 6* ^ 0. 



In the following we shall use the notation k = 3 In 



this is a slow-varying. 
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empirical quantity, which varies only slightly as the droplet spreads and introduces 
minor, logarithmic corrections to the power laws; experimental data suggest that a 
good choice is a nearly-constant k « 120 [ISj. 

We now come to the dissipation in the film. A striking result of Hervet and de 
Gennes [11] is that the complete wetting regime is characterized by 

Wr^^so^SU, (4) 

which means that the dissipation within the film compensates exactly the first term 
on the r.h.s. of rendering the rate of spreading independent of S. This result was 
obtained for non-retarded van der Waals substrate forces, but can be generalized. 

Finally, the form of the dissipation term Wmicro was discussed by Blake and 
Haynes [291, '30]: it was found that Wmicro ~ CU"^ at leading order in U, where C is a 
constant friction coefficient. Note that in the case of complete wetting C is dependent 
on the thickness of the precursor. 

Thus, provided that S is consumed entirely in the film, the dynamical behavior 
results from a competition between the two remaining dissipation channels, Wmacro 
and Wmicro- As — > 0, Wmicro IS ^-independent, whereas Wmacro ^ 1/6* and thus 
clearly dominates at long spreading times. Dissipation at the microscopic contact line 
may dominate (e.g., for low viscosity liquids) at intermediate times, but ultimately 
hydrodynamic dissipation in the core drop will take over [U [H] . 

Consequently, neglecting Wmicro as compared to Wmacro, one finds in the small-0 
limit that HI) adopts the following form: 

6l^w2KCa, (5) 

which is a fundamental relation between the velocity of the moving contact line and 
the instantaneous value of the contact angle (Ca — r]\J/<7 is known as the capillary 
number) . It was derived analytically by Voinov [4] and by Tanner [5] using a different 
approach (within the lubrication approximation). 

Now taking into account that the volume V ~ jR^9 of the bulk drop remains 
nearly constant during spreading, we obtain 

from which the Tanner laws R ~ t^/^^ and 6 ^ t^'^/^" ensue trivially. The behavior 
described by ^ and iJH) has been observed experimentally in [5| fT2lfT3] . 



3. Energy balance approach revisited 

Equations |(2]), |(5]) and ([6]), under the assumption of well separated length scales of 
the bulk and film, predict an unbounded growth of R. This is an ideal spreading 
behaviour, through which the droplet virtually thins down to a molecular film. As we 
have already remarked, this is not always the the spreading may terminate 

with the appearance of equilibrium pancake-like structures [U [6l [7]. In the latter 
situation, the Tanner law clearly describes an intermediate stage, and the transition 
to a pancake must be described by some kind of crossover in terms of R{t): our 
intuition is that R will tend to a finite value although it is not clear whether the 
definition of R will remain consistent with figure la. 

An opposite trend was revealed by recent experimental studies focusing on the 
spontaneous spreading of nematic liquid crystals (cyanobiphenyl 5CB) on hydrophilic 
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[9] or hydrophobic [l^ substrates: after a transient Tanner stage, an "acceleration" of 
the spreading process has been observed. The data suggest that the base radius R, as 
inferred from the inflection point of the thickness profile, grows algebraically but with 
an exponent which is substantially larger than — 0.1. In ^ it was shown that the 
Tanner law crosses over to R ^ t" with a ~ 0.2. Later it was realized (see figure 6 
in [T^) that the Tanner relation in ^ does not hold for late stages of spreading: for 
small 9 and Ca the best fit to the experimental data follows 9 '-^ Ca°' rather than that 
predicted by The latter relation, together with the volume conservation condition 
R^9 ~ V yields R ^ t"' with a ~ 0.19. A more thorough analysis of the behavior 
depicted in figure 5 in [10] suggests that actually the reported law R ~ t^^^ is only a 
part of a crossover from the Tanner stage to an even faster growth law. Experimental 
data in [T^ span time scales ranging from a second to two hours, and at the end of 
the experiment the trend is clearly rather R ~ t^/^ than R ~ i^/^, and possibly still 
accelerating. Consequently, although there is no conclusive evidence on the precise 
value of the exponent a characterizing the accelerated spreading regime, it is clear that 
a is significantly larger than the Tanner exponent and thus the physical mechanism 
responsible for the late, post-Tanner stages of spreading might be different from the 
one described in section 2. 

3.1. First guess: shear thinning 

We notice first that the accelerated spreading was observed for nematic liquid crystals. 
Nematic crystals are known to have a non-Newtonian, shear-thinning rheology (see, 
e.g., [HH [32]). Shear thinning affects the flow pattern, which necessarily modifles 
the spreading dynamics. Thus our flrst idea is to revisit the l.h.s. of HI) and, more 
specifically, the term Wmacro in The expression ([4]) for Wmoso is also queried in 
the case of shear thinning. 

A detailed analysis of the contact line dynamics within the framework of the 
thin film model shows that the characteristic shear rates in the capillary wedge and 
in the film decrease as the contact line velocity decreases p3]. In consequence, for a 
spontaneously spreading droplet of a non-Newtonian, shear-thinning fiuid the effective 
viscosity will increase with time resulting in a spreading law of the form R ^ t"' 
with a < 1/10. Numerical simulations carried out in [34] confirm that a < 1/10 
for shear-thinning fiuids and that a > 1/10 for shear-thickening fiuids. Hence the 
dominant effect from shear thinning is that the spontaneous spreading of a non- 
Newtonian fiuid is generally slower than predicted by the Tanner law and can not 
explain the experimentally observed acceleration of the spreading process. 

3.2. Second guess: line tension 

We now turn our attention to the r.h.s. of ([T]) and notice that the unbalanced Young 
force - which is also the r.h.s. of ([2]) - is, in fact, a mere approximation of the actual 
driving force of spreading. In general, the total free energy of the Hquid/soHd system 
can be decomposed into bulk, surface, line, and point contributions (see, e.g., |17[ll8j). 
Thus, the driving force ([2]), as a derivative of the total free energy, should also contain 
all these contributions. This picture, of course, is meaningful only in the case of 
macroscopic drops, for which both the surface tension and the line tension are well 
defined. 

As a matter of fact, the spherical cap adequately describes the profile of 
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the bulk drop, but the profile of the mesoscopic film deviates from it, such that 
the quasistationary free energy 2nRF must include a correction term, which is 
accumulated in the vicinity of the apparent contact line and can be seen as a line 
energy r multiplied by the apparent perimeter 27ri? [35]. The expression ^ takes 
into account the surface energies of the three macroscopic interfaces meeting at the 
apparent contact line, but does not include the line tension contribution. 

The idea that the Hue tension may have an appreciable impact on the global 
behavior is not new. As an excess quantity, r can be positive or negative, as noticed 
already by Gibbs [36]. Negative line tension, for example, can significantly reduce the 
work required to create a nucleus (100 Angstrom in diameter) of a new phase on solid 
or liquid substrates [37]. Conversely, positive line tension can explain the stability of 
Newtonian black films towards rupture [38]. For liquid droplets of nanometer size, 
negative (resp. positive) fine tension can promote spreading (resp. dewetting) even 
if the macroscopic spreading parameter S is negative (resp. positive) [35]. A review 
of different phenomena caused by line tension effects and some conceptual aspects of 
line tension can be found in |17[ 118). 

Evaluation of the contribution due to the fine tension r involves many delicate 
issues (e.g., a proper definition of the effective interface potential used in the model, or 
a proper convention when choosing the Gibbs dividing interface) and, in general, is a 
more complex problem than the calculation of the surface tension - essentially because 
more phases meet at the contact line than at an interface (see, e.g., \T7\ [TS] for a 
more thorough discussion). The problem is already difficult in equihbrium situations 
(e.g., partial wetting, with S < 0), and clearly becomes even more complex when 
one considers spontaneous spreading, since here one has to account for the temporal 
evolution of the droplet thickness profiles. 

The consideration of these subtle points is beyond the scope of the present 
approach. For our purposes it will be sufficient to resort to a recently proposed 
phenomenological generalization of ((H) in terms of non-equilibrium thermodynamics: 
as shown in |39], the force applied to the apparent contact fine of a droplet can be 
written down as 

fr = S + <j{l-COS0)-^, (7) 

which differs from the expression in ^ by an additional term accounting for the 
contribution of the contact line tension r to the driving force of spreading. The 
definition of fr as a generalized Young force is valid both in complete and partial 
wetting, and is consistent with the so-called modified Young equation fr = 0, obtained 
at equilibrium by an appropriate generalization of Gibbs classical theory of capillarity 

IMIIID]. 

At this point we must stress that the expression ^ for fr is formally valid only 
if T is constant. If we assume a power-law behaviour for r ^ , then the modified 
Young force becomes S + a {1 ~ cos 9) — {1 + P). Here we argue that it is not likely 
for T to vanish at long t and large R, and thus f3 > 0. Thereby ^ essentially holds 
for r ~ R^, up to a numerical factor {1+ P) > appHed to the fine tension term, i.e., 
the last term on the r.h.s. of ([7]). This consideration has little impact on the following 
qualitative argument, but will be relevant to quantitative implications. 
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3.3. Line tension effects on spreading 

Suppose now that HI) holds; that the dissipation Wmoso in the film obeys ([4]); that 
the dissipation Wmacro is given by ^ ; but that the driving force of spreading is now 
determined by |[7|, i.e., that the line tension contribution is taken into account . Then 
((5|) is replaced with the following: 

Since t is typically very small, one naturally finds that the surface tension 
contribution will dominate at small and intermediate times, again giving rise to the 
Tanner law R ^ t^/^°. On the other hand, as R grows to sufficiently high values, 
the line tension contribution will inevitably take over and become a dominant driving 
force, provided that t/R decays slower than the capillary term aO^. 

In the latter regime, the behavior is crucially dependent on the sign of t: 

• H T is positive and tends to a constant value, which is physically plausible, then 
llH) predicts that spreading will terminate at a finite value of R, for which the first 
and the second terms on the r.h.s. of ^ become equal to each other. One may 
interpret this as an indication of the formation of a pancake. However, during the 
actual transition to a pancake, a drop would not retain the shape of a capillary 
cap (see figure lb), which somewhat challenges this prediction. 

• If r is negative and the second term in ^ dominates, we find the following 
post-Tanner behavior: 

R^ (~ f Tdt] . (9) 



This spreading law is qualitatively different from R ~ t^/^° . In the framework 
outlined in section 2, the driving force of spreading is associated with the surface 
tension. By contrast, during the late stages of spreading, the droplet becomes 
flatter and can be viewed as "quasi two-dimensional". It is then not surprising that 
the Hne tension r should govern the spreading process, provided that |t| 3> tjRO"^. 

It is tempting to obtain coarse quantitative results from ([9]) and from the condition 
|t| ^ (jR6^ . In particular, if we assume that r is a negative constant, then ([9]) predicts 
i?~ which agrees with previously reported experimental results [9l[l0]. Looking 
at figure 5 in [TO] we can estimate the value of r from the characteristic base radius 
at the apparent crossover between the Tanner stage and the accelerated spreading 
regime: this yields r « — 10~^ N. This value is an order of magnitude higher than 
previously reported values of the line tension in the partial wetting situations, but 
it must be noted that many experimental measurements of r have been performed 
for simple Hquids; in the case of nematic liquid crystals an elastic contribution to 
the effective interface potential (a consequence of the anchoring properties) may yield 
substantially higher values of r. 

However, upon a closer examination of the case of negative t, there is no good 
reason to expect that r should approach a constant value. One rather expects that 
the film region will become progressively more pronounced and the drop profile will 
significantly deviate from a spherical cap-like shape; hence r, a functional of local 
droplet thickness, will grow as a function of time (in terms of its absolute value). 
Consequently, at late spreading stages, one may expect a growth of R{t) that is faster 
than R--t^/^. 
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4. Conclusions 

We have presented in this paper both the classical energy balance approach - as 
developed by de Gennes - and a revised version of it, which incorporates a line tension 
contribution to the driving force of spreading. The revisited framework was motivated 
by apparently contradictory trends at long spreading times for macroscopic droplets 
in complete wetting. By taking line tension into account, we have complemented the 
classical framework with the following twofold interpretation: 

• A positive line tension - essentially a "collar" around a spreading droplet - stops 
spreading and is responsible for the formation of mesoscopic pancakes. 

• A negative line tension - which tends to lengthen the apparent contact line - 
governs the late stages of spreading, resulting in the acceleration of this process. 

We must now voice a few words of caution concerning our approach, which seems 
intuitive but may have several shortcomings due to its simplicity. 

First we admit that the formal definition of line tension t in partial wetting - a 
functional of an equilibrium profile, and an integral of the effective interface potential 
- cannot be easily generalized to complete wetting and thus remained undetermined 
within our analysis. The case of constant, positive r is plausible, but negative r is 
more likely to grow as a function of time and base radius, a growth which we are 
unable to specify. We merely expect that, as in the partial wetting case, r is a certain 
functional of the mesoscopic thickness profile, rather than an independent, arbitrary 
quantity (see, e.g., |41|). 

We have already stated that in the case of positive r our prediction of the 
emergence of pancakes is indicative at best: indeed, in the process of reaching the 
stationary shape of a pancake, a droplet gradually deviates from the spherical cap 
shape assumed by sections 2 and 3; this deviation entails corrections that our approach 
does not account for. A similar word of caution exists for negative r and arises from a 
thorough comparison with experiment. Indeed, a remarkable feature of nematic 5CB 
droplets observed in [D', ToJ is that the reported accelerating phase is accompanied by 
the development of a large "foot" (essentially the drop adopts a bell shape similar to 
figure la, without exaggeration). This large foot is a warning against the applicability 
of energy balance as developed in sections 2 and 3: the key assumption of a spherical 
cap of constant volume Vcap = jR^0 is less and less valid as the mesoscopic film drains 
liquid from the macroscopic droplet. In other terms, the separation of macroscopic 
and mesoscopic length scales (both vertical and lateral), in the experimental layout 
that we are trying to describe, may be more precarious than what was assumed in our 
analytical framework. 

In the light of these shortcomings, our agenda is to develop a more robust 
approach describing late stages of droplet spreading, based on the seminal approach by 
Tanner describing the time evolution of the thickness profile of whole droplets within 
the thin film approximation [H [5] . In the latter framework, scale separation is not an 
issue, and the fundamental result is that droplets gradually turn into diffusive films in 
the sense of Derjaguin [42]. Results have already been obtained for the specific case 
of nematic droplets, and we shall present them in a companion paper j43j . 

We are also looking forward to developing the notion of a dynamic line tension in 
spreading processes, through a detailed study of the hydrodynamic wedge as in [T|[ii]. 
It is important to note that nematic 5CB droplets exhibit both of the non-trivial types 
of post-Tanner behaviour detailed in this paper: R{t) crosses over from R ~ t^^^ to 
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faster power laws (at this point the drop looks like figure la), but eventually the 
spreading terminates with a mesoscopic pancake (figure lb). In terms of r, this means 
that the physically relevant line tension is negative during the accelerating phase, but 
later becomes positive. As figure 1 suggests, the geometrical properties of the droplet 
are quite different in both regimes, and it is not clear whether a consistent definition 
of T can be established for such liquids as nematic Hquid crystals. 
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